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Low-temperature transport experiments on a p-type GaAs quantum dot capacitively coupled to 
a quantum point contact are presented. The time- averaged as weU as time-resolved detection of 
charging events of the dot are demonstrated and they are used to extract the tunnelling rates into 
and out of the quantum dot. The extracted rates exhibit a super- linear enhancement with the 
bias applied across the dot which is interpreted in terms of a dense spectrum of excited states 
contributing to the transport, characteristic for heavy hole systems. The full counting statistics of 
charge transfer events and the effect of back action is studied. The normal cumulants as well as the 
recently proposed factorial cumulants are calculated and discussed in view of their importance for 
interacting systems. 
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I. INTRODUCTION 

Quantum dots (QDs or simply dots) are small con- 
ducting islands containing a tunable number of confined 
charge carriers tunnel-coupled to source and drain leads. ^ 
This tunability is achieved by means of an additional 
capacitive coupling of the QD to nearby gates. It is a 
fascinating experimental observation that a similar ca- 
pacitive coupling to a nearby electrical current passing 
through a constriction provides the possibility of mea- 
suring the charge of the dot with a precision of a small 
fraction of an electron's charge.^ The conductance of the 
constriction changes as a function of the average charge 
population of the QD. 

Two-level fluctuations (random telegraph noise) of the 
detector current provide more information about the 
dot than just the average current. The fluctuations of 
the current enable the time-resolved detection of single- 
particle charging and de-charging events in the QD.^ ^ 
This, on the other hand, reveals more information about 
the energy spectrum of the dot, the relaxation of ex- 
cited states to the ground state and their coupling to the 
leads. ^ This information is valuable for the case of hole 
QDs where the present understanding is limited by the 
lack of experimental results. 

Counting statistics of the charge transfer is another 
tool to study quantum dots. Experimental studies of 
counting statistics using charge detection with a quan- 
tum ^int contact (QPC) were started by Gustavsson 
et al.,^^ and Fujisawa et al.,^ and continued by Fricke 
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et al.12^ All these experiments were performed on n- 
type GaAs or InAs electronic systems.^ Therefore it is 
interesting to compare these results with those obtained 
in a QD realized on a p-GaAs two dimensional hole gas 
(2DHG), where the carrier-carrier interactions are sup- 
posedly stronger both in the dot and in the leads com- 
pared to n-type systems. 

In this article we investigate these effects in a p-type 
GaAs QD system. The large effective mass of holes, sev- 
eral times larger than that of conduction band electrons 
makes carrier-carrier interaction effects more pronounced 
compared to the kinetic energy than in their electronic 
count erparts.fl^ The same reason leads to the fact that 
screening is expected to be stronger and that the single- 
particle energy spacing is much smaller in confined p-type 
systems, making it very difficult to be resolved at acces- 
sible temperatures.^^ Additionally, the strong spin-orbit 
interaction in the valence band holds promise for inter- 
esting spin physics in these QDs.^ 

However, the fabrication of tunable p-type QDs is chal- 
lenging,^^ essentially because metallic gates on top of 
shallow p-doped heterostructures have a low Schottky 
barrier resulting in leaky and hysteretic behavior and 
different fabrication techniques have to be adopted. The 
Coulomb blockade effect in lithographically-defined dots 
in p-type GaAs heterostructures was first demonstrated 
by Grbic et using local oxidation lithography. The 
same technique was later shown to be effective in fur- 
ther confining the carriers and observing the individ- 
ual excited energy states in a QD.I^ Induced SETs were 
also fabricated using undoped GaAs heterostructures and 
were shown to exhibit Coulomb blockade effects. In 
spite of this progress the level of control on the fabri- 
cation of these nanostructures and the understanding of 
the role of interactions in hole systems is still far from 
complete. In this article we attempt to improve on this 
understanding by realizing time-resolved charge detec- 
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tion of hole tunnelling into a QD fabricated by shallow 
wet chemical etching. 



II. SAMPLE AND SETUP 

Fig.jlja) shows an AFM micrograph of the sample 
which was patterned in the 2DHG by electron beam 
lithography followed by shallow wet chemical etching!^ 
The trenches seen in Fig.jlja) are 20 nm deep and lo- 
cally deplete the 2DHG situated 45 nm below the surface, 
thereby separating the 2DHG plane into laterally discon- 
nected regions. Each of them is connected to metallic 
leads via ohmic contacts. The host material consists of 
a C-doped GaAs/AlGaAs heterostructure grown along 
the (100) plane. Prior to sample fabrication the qual- 
ity of the 2DHG was characterized by standard magne- 
totransport measurements at 4.2 K and a hole density 
of n = 2.7x 10^^ cm-2, and a mobility of /i = 60'000 
cm^/Vs were obtained. 

The sample consists of a QD together with a nearby 
QPC. The measurement setup and the applied bias volt- 
ages are also schematically shown in Fig.jlja). The dot 
bias and the QPC bias are both applied symmetrically. 
The overall potential of the QPC (Vqpc) is used to con- 
trol the electrochemical potential of the QD, while the 
plunger gate (PG) is used to tune the QPC transmis- 
sion. The in-plane gates Gl and G2 are used to tune 
the tunnel coupling between the QD and source (S) and 
drain (D), but they also have a significant lever arm on 
the dot. 



III. RESULTS AND DISCUSSION 
A. Time- averaged/ time-resolved charge detection 

Fig.jljb) shows simultaneous measurements of QPC 
and QD currents as a function of the voltage applied to 
the gate G2 at the temperature T ^ 1.2 K. As the gate 
voltage is increased, the holes are unloaded from the dot 
one by one. The dot current shows clear conductance res- 
onances at the charge degeneracy points where the charge 
state of the dot changes by one elementary charge. This 
can be clearly seen as a step of 30 pA in the QPC current 
4%) at the position of the Coulomb peaks. Note that 
the average QPC current decreases with G2 (due to the 
corresponding lever arm), since no electrostatic compen- 
sation was performed here in order to avoid activating 
additional fiuctuators in the sample which degrade the 
detector signal. For the remainder of the paper we will 
consider the results obtained in a dilution refrigerator 
with a base temperature of 100 mK. 

When the bandwidth Tjj of the detector circuit is small 
compared to the tunnelling rates of the QD, it only re- 
sponds to the average charge population of the dot. As 
F/) is increased to about 3 kHz and the tunnelling bar- 
riers are tuned sufficiently opaque, time-resolved charge 
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FIG. 1: (a) AFM micrograph of the sample consisting 
of a QD with an integrated charge read-out QPC. The dark 
regions are trenches created on the surface of the heterostruc- 
ture by chemical etching. The applied voltages are shown 
on the same figure, (b) QD current (blue) and QPC current 
(green) as a function of Vg2 measured at the temperature of 
T ^ 1.2 K. (c) QPC current as a function of time, showing a 
few holes tunnelling into and out of the QD over a timescale 
of 200 /is. The lower current level corresponds to a state when 
the dot holds one excess hole. The QPC current was filtered 
with a 3 kHz software filter and re-sampled at a frequency of 
14 kHz. The random variables Tin and Tout quantify the times 
it takes for a hole to tunnel into and out of the dot, respec- 
tively, (d) The normalized histogram of the detector current 
showing two distinct current levels corresponding to the two 
charge states of the QD (e) The histogram of times a hole 
needs to tunnel into the QD (blue), tunnel out of the QD 
(green), and the event time defined as the sum of two con- 
secutive tunnelling in and out events (red) for a symmetric 
configuration (Fin = Tout)- The green curve exactly overlaps 
with the blue curve, (f) The occupation of the QD and the 
number of charge events as the gate voltage is swept over a 
Coulomb blockade peak with fits to the Fermi-Dirac distribu- 
tion (green curve) and its derivative (red curve), respectively. 



detection becomes possible. The detector current then 
exhibits a two-level ffuctuating behavior as a function 
of time because holes tunnel into and out of the dot 
[shown in Fig.jljc)]. The two levels on the histogram 
of the detector current [Fig.jljd)] are a result of the 
two charging states of the dot. The random variables 
'?"in/out quantify the time it takes for a hole to tunnel 
into or out of the QD and are used to calculate the 
tunnelling rates according to Fin/out = ('^in/out )~"'^ (Here 
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bracket denotes an ensemble averaging). The latter can 
be used to quantify the coupling symmetry of the QD 
to the leads by defining the normalized coupling asym- 
metry a = {Tin - rout)/(rin + Tout). The histograms of 
^in/out and the event length (revent = ^in + Tout) are plot- 
ted in Fig.jlje) for a symmetric configuration of the QD 
(a = 0). About 2 million events, accumulated over more 
than 5 hours, were used to produce these histograms in- 
dicating the stability of the sample. The exponential 
distribution of the tunnelling rates motivates the use of 
the rate equation technique to study the statistics of the 
charge transfer. The short-time suppression of Tevent is 
a consequence of the correlated transport and sequential 
tunnelling through the QD.^ 

Figure [TJf) shows the occupation probability of the 
QD (extracted from the duty cycle of the detector cur- 
rent) together with the average number of events (in 1 
second time traces) as the gate voltage is swept over 
a Coulomb blockade peak. They can be fitted with a 
Fermi-Dirac distribution, from which a hole temperature 
of Thole ^ 180 mK is obtained. 



B. Excited state spectrum 

Fig.|2|a) shows the event rate as a function of the QD 
bias VsB and the gate voltage Vg2' Due to a strong 
energy-dependence of the tunnelling rates, tunnelling on 
the adjacent Coulomb peaks is either too fast or too slow 
to be properly detected. The charging energy of the QD 
is Ec ^ 2 meV which corresponds to a total capaci- 
tance of Cy: ~ 80 aF. Assuming a disk-like shape for the 
dot with Cy: = SeoErr^ where r is the dot radius and 
Er = 12.9 for GaAs, this provides an (upper) estimate 
of ~ 170 nm for the electronic diameter of the dot, and 
an upper limit of 90 for the number of holes in the QD. 
With this diameter the mean single-particle level spac- 
ing can be calculated from A = Tih^ /m*A with A = 7rr^, 
giving A ^ 20 /iV comparable to ksT. The large ef- 
fective mass of the holes results in a dense spectrum of 
confined states and therefore the hole QD can be consid- 
ered to be in the crossover between electron QDs with 
a discrete, and metallic SETs with a continuous excited 
state spectrum for the sizes investigated here. This man- 
ifests itself in the fact that it is not possible to resolve 
excited states in the diamond measurements as shown in 
Fig.|2j Were this resolution possible, we would expect a 
stepwise increase of the number of events with the steps 
parallel to the edges of the diamond.^ Nevertheless it can 
already be seen in this figure that the number of events 
generally increases with increasing QD bias. Note that 
the event rate in Fig.|2|a) is not symmetric with applied 
bias, presumably due to some tunnelling asymmetry of 
the barriers. 

More insight into the role of the dot excitation spec- 
trum is obtained by looking at Fig.[2|b)-(c) where Fin 
and Tout are plotted instead of the number of events. 
The increase in Fin and Fout with Vsd can be understood 
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FIG. 2: (a) Event rate Fevent, and tunnelling rates (b) Fin 
(c) Fout as a function of the gate voltage and the applied QD 
bias demonstrating half of a Coulomb diamond. All the rates 
increase with the dot bias. White arrows indicate the relative 
position of the QD electrochemical potential ii% with respect 
to those of source and drain, (d) In the case of a moderate 
energy-dependence of barriers, tunnelling into and out of the 
dot both involves many excited states and relaxation does not 
change the qualitative picture. 



in terms of the additional available tunnelling channels 
in a system with dense spectrum. The relative position 
of the QD electrochemical potential /i^^ with respect to 
those of left {(11) and right {(Ir) leads are indicated with 
white arrows. Lines of constant Fin/out are parallel to 
the edges of the diamonds. In particular. Fin depends 
only on the difference of the electrochemical potentials 
of the source {(is) and the dot (is — Mat (^^^ positive bias 
IJ^s = l-^L and /J^D = I^R while for negative bias /i^ = /j^r 
and = jiiL)- This suggests that the number of avail- 
able (excited) states between these two levels is the cause 
of the increase in the tunnelling rate. Provided that they 
have enough energy, the holes tunnelling-in can occupy 
any of these states and this increases the tunnelling rate 
with bias. Similarly, Fout depends only on the difference 
between the electrochemical potentials of the drain {/j^d) 
and the dot /i^^ — I^d^ meaning that the number of op- 
tions for holes tunnelling-out also increases with the bias. 
For example, it is possible for a hole in the {N + l)-hole 
ground state to tunnel out and leave the dot in any of 
the A/'-hole excited states. Motivated by these ideas, the 
level diagram of the dot is represented in Fig.|2|d) with 
a dense ladder of excited states both above and below 
the ground state transition /i^^ with an electrochemical 
potential of /i^^ (the index m refers to a transition in- 
volving excited states) which contribute to Fin and Fout, 
respectively. 

An alternative explanation for the bias-dependence of 
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Tout would be to assume a strong energy-dependence of 
the tunnelling rates together with a ladder of excited 
state transitions only above the ground-state transition 
and a slow relaxation rate. Since the dot has an almost 
continuous spectrum the holes can tunnel in and out elas- 
tically at all energies. In the presence of a strong energy- 
dependence of the tunnelling rates, the higher levels can 
have much faster tunnel rates. Therefore tunnelling-in 
happens basically within a ksT window around fis (op- 
timized by the corresponding Fermi distribution) . Due to 
slow energy relaxation in the dot, tunnelling-out would 
also happen at the same energy and therefore increase 
very fast with the bias in agreement with the data. More- 
over, the tunnelled-in holes would have a finite chance to 
tunnel back to the source. The measurements of the full 
counting statistics that are provided in the next sections, 
however, point to the presence of a relaxation rate much 
higher than the bandwidth and therefore excited state 
transitions only above the ground state transition can- 
not explain the measurement results. 



C. Rate equation simulation 

To verify the ideas discussed in the previous section we 
did rate equation calculations, first for a dot with a sin- 
gle state in the A/'-charge configuration and few excited 
states in the (A^ + l)-charge configuration and second, for 
the case where both N and (A^ + l)-charge configurations 
have many excited states. The occupation probabilities 
{pf^ and pf~^^) of individual states are calculated in the 
steady state and the total tunnelling-in/out rates are ob- 
tained from (no summation on i for the first model) 



J- out — J- N^N+l/^j / / ,P 
ij 3 

ij i 

The parameters of the model (/(e) = 1 — /(e)) 



(1) 



^ N+l^N 



depend on the gate and the applied bias through the ar- 
gument of Fermi distributions /(e). All the states in the 
dot are assumed to be coupled to the leads with the same 
coupling (Ti^ =rii = 100 Hz). Since we expect our QD 
to be far from the few- hole regime, a linear spectrum with 
a constant mean-level spacing is assumed: 20 levels with 
the energy- separation of 100 /ieV are taken into account. 
A strong relaxation {R = 1 kHz) is assumed for all levels. 

In the first model in which there are no excited states 
for the A^-charge configuration. Fin increases each time 
an excited-state transition of the (A^ + l)-charge configu- 
ration enters the bias window. However, Fout stays con- 
stant and while the event rate also increases with bias, it 
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FIG. 3: Rate equation simulation of a dot with a dense 
spectrum of 100 /xeV equi-distant energy levels. Fin in (a) 
and Fout in (b) both increase with the dot bias. The equi-Fin 
lines and equi-Fout lines are parallel to source and drain lines, 
xrespectively. (c) Fevent- (d) Horizontal cut through Fin, Fout 
and Fevent in figures (a)-(c) at resonance, (e) Same as (d) but 
assuming an exponential energy-dependence of barriers, (f) 
Horizontal cut through the measurement data in Fig.[2|a)-(c) 
at resonance as in (d). 



saturates at the value of the tunnelling-out rate, which 
becomes the bottleneck. 

In the second model, motivated by the discussion in 
the previous section, a similar spectrum is assumed for 
the A/'-charge configuration. This time Fout also increases 
with the bias and this simple model is able to qual- 
itatively reproduce the measurement result, as shown 
in Fig. 3 Fig. 3rd) shows a horizontal cut through fig- 
ures|3|ay-(c) at resonance. The linear increase of the tun- 
nelling rates with bias is a consequence of equal tunnel 
couplings of individual states which is slightly different 
than the non-linear increase of the tunnelling rates ob- 
served in the measurement (Fig.jsjf)), but it must be 
noted that the difference can be easily captured by as- 
suming an energy-dependence of the tunnelling rates as 
shown in Fig.jsj^e). 



D. Counting statistics 

Fig.Qa) shows the histogram of the number of events 
in a 50 msec time window with symmetric tunnelling cou- 
pling (a = 0) of the dot (point A in Fig.[2|a)). Three dis- 
tributions, namely the Poisson distribution, the Gaussian 
distribution and the model of Bagrets-Nazaro\^^ with 
two-states are plotted in the same figure. The fact that 
the model of Bagrets-Nazarov (BN) matches perfectly to 
the data indicates that, in spite of the dense spectrum 
of the dot and its contribution to transport (as shown 
in the previous section), the statistics is dominated by a 
two-state Markovian model. This is presumably due to 
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FIG. 4: (a) Histogram of the number of events during 
T=50 msec at the point A in Fig.[2|a) with symmetric barri- 
ers (a = 0). The blue and green curves show the Poisson and 
Gaussian distributions respectively, calculated using the mean 
and variance of the measured data (all these distributions are 
discrete and the connecting lines are just guides to the eye). 
While the Gaussian distribution fits reasonably well to the 
histogram, the data is best described by the Bagrets-Nazarov 
distribution with the input parameters Fin and Fout- (b) The 
first two normalized cumulants of the statistical distribution 
of events as a function of tunnelling barrier asymmetry a, 
calculated from 100 time traces, each 10 sec long. The blue 
points show C2/C1 or the Fano factor and the red points 
show C3/C1 which is the skewness. The solid lines are the 
model predictions. While the Fano factor agrees quite well 
with the model the skewness points are scattered much more 
due to limited statistics, (c) Fin, Fout and Fevent as a func- 
tion of the bias voltage on the QPC. The decrease m i in and 
increase in Fout is most probably due to a gating effect, (d) 
Fano factor and skewness as a function of QPC bias showing 
that the statistics of electron transport is not influenced by 
the emission of energy quanta by the QPC. The dashed lines 
show the model predictions discussed in the text augmented 
by finite-bandwidth corrections (B.W.C).'^ The red shaded 
area shows the onset of the detector signal degradation due 
to charge fluctuators in the QPC (Counting was not possible 
for Fc^jfg < 50 fiV). 



a strong relaxation in the quantum dot. 

The Fano factor {F = C2/C1) and the skewness {S — 
C3/C1) extracted from the data are plotted as a func- 
tion of the asymmetry a in Fig.jljb) together with the 
corresponding predictions of the BN model {Cn is n-th 
cumulant of number of events in a series of 10 sec- long 
time traces). The agreement of the model with the data 
indicates that again a two-state Markovian model is suf- 
ficient to describe the observed statistics. 

It is not a priori clear whether each event in a given 
time trace corresponds to a charge transfer from source 
to the drain. In a quantum dot without any excited 



states in the relevant energy window, the ratio between 
the charges tunnelled back to the source to those trans- 
ferred to the drain is essentially ksT/YsD which is 1/40 
at the point A in Fig.[2|a). However, in a quantum dot 
with a dense spectrum with energy-dependent tunnelling 
rates the ratio can be higher. In the presence of a strong 
energy-dependence of the tunnelling rates, relaxation is 
crucial to ensure fully uni-directional transport. While 
this condition Fin/out <C 7 is presumably satisfied in 
our case (Fin/out Fi) <C 7) with the bandwidth of 
T]j = 3 kHz, this might not be the case in more strongly 
coupled regimes. 



E. QPC back-action 

The power dissipated in and around the QPC is emit- 
ted as photons and phonons close to the QPC and hence 
may cause back-action on the QD either by increasing 
the effective temperature of the leads, or by excita- 
tion of the QD due to photon and phonon-assisted tun- 
nelling (PAT).^^ These PAT effects are usually under- 
stood in terms of energy transfer between the QPC and 
the electrons/holes in the dot so that they could over- 
come the relevant energy barrier (Coulomb blockade or 
single-particle level spacing). Therefore they are char- 
acterized by an energy cut-off corresponding to the mean 
level spacing of the dot, below which this energy transfer 
does not take place. Identifying the dense spectrum of 
our hole QD as the source of the peculiar bias depen- 
dence of the tunnelling rates, it would be interesting to 
see if the detector has any back-action on the dot due to 
PAT and how much it contributes to the transport and 
its statistical properties. 

Fig.[4](c) shows how the variation of the bias on the de- 
tector QPC influences the tunnelling rates of the dot. For 
small QPC bias (Vqpc < 70 /J.V) detection is not possible 
due to low signal-to- noise ratio while for large QPC bias 
(^QPC > ^00 /iV) many fluctuators in the QPC are acti- 
vated and the overall quality of the signal is degraded by 
the additional telegraph noise due to these fluctuators. 
The red shaded area shows the onset of this degradation. 
Fig.[4]^d) shows the effect of QPC bias on the Fano fac- 
tor and skewness of the hole transfer distribution. The 
agreement between the measurements and the two-state 
Markovian BN model (dashed line) implies that the ef- 
fect of the QPC on the dot can be phenomenologically 
lumped into the tunnelling rates Fin and Fout- This is 
in contrast to what is expected from a master equation 
calculation, which can be slightly modified to include the 
effect of PATs. It is shown in the Appendix using a simple 
model that the presence of an excited state and photo- 
excitation of a QD generally alters the statistics obeyed 
by a dot unless relaxation is faster than both the tun- 
nelling rates and the excitation rate. 

Furthermore, while Fout shown in Fig.|4jc) increases 
monotonically with the QPC bias. Fin decreases, sug- 
gesting that the influence of QPC on the dot is at least 
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partially a simple gating effect. Considering the close 
proximity of the QPC leads and the dot tunnel barriers 
in Fig.jlja) and the equal polarity of the dot and QPC 
biases, this is not surprising as most of the applied bias 
voltage drops over the QPC. As a result, the height of 
the source tunnelling barrier increases (decreasing Fin) 
and the height of the drain tunnelling barrier decreases 
(increasing Fout) for positive QPC bias. Moreover, since 
a symmetric bias of 700 /iV is applied to the dot (point A 
in Fig.[2|a)), Fout is expected to exhibit a step at a QPC 
bias of about 350 /iV as the holes can tunnel out to the 
source lead, while the measurement shows a monotonic 
increase of Fout • Similar experiments were performed in 
an off-resonance configuration and no change in the tun- 
nelling rates were observed for Vqpq < 600 /iV. Therefore 
we conclude that the relaxation is dominant in our exper- 
iment and the results of Fig. 4]^c)-(d) are mainly gating 
effects in the window of QPC biases investigated. For the 
other counting measurements the QPC bias was kept at 
250 /iV. 



struction. Recently, Kambly et al.'^ proposed the use of 
factorial cumulants (FCs) for this purpose as any zero- 
crossing oscillations in the latter directly indicate the off- 
set of ZGF from the real axis pointing towards relevance 
of interactions. We have calculated the first twelve FCs 
from our data which are shown in Fig.|5|c)-(d) and on the 
logarithmic scale in Fig.|5|e)-(f). Due to the logarithmic 
scaling of the FCs the latter plot is more convenient to 
follow the evolution of the results. Note that consec- 
utive FCs alternate sign as indicated by red and blue 
colors. For a two-state Markovian system no oscillations 
in the factorial cumulants are expected in agreement with 
the fact that there is no clear zero-crossing oscillations in 
the data. This again implies that the two-level system 
is a surprisingly good model to describe the statistical 
properties of our multi-level QD presumably due to the 
strong relaxation, as explained in the Appendix. 



IV. CONCLUSION 



F. Normal vs. factorial cumulants 

For a closer look at the statistics, we have calculated 
the first twelve cumulants of the tunnelling events and 
plotted them together with the predictions of the two- 
state Markovian BN model in Fig.[5|a)-(b). In gen- 
eral the finite bandwidth, the limited signal-to-noise ra- 
tio and the finite statistics influence the calculation of 
the cumulants. While the first two problems can be in 
principle taken into account by introducing additional 
Markovian states into the model, the finite statis- 
tics is responsible for the error bars in Fig.[5|a)-(b). 
The latter is calculated from the covariance formula^ 
{ACnACm) = mla^'^SmnN-^ + 0(7V-2). The N in the 
denominator signifies the importance of the amount of 
statistics for a reasonable accuracy. For a fixed total 
number of events K (two million in our case) used to 
calculate the cumulants as a function of (n) , the amount 
of statistics is equal to N = K/ (n). Also note that 
C2 = cr^ eventually grows linearly with (n) in the steady- 
state and therefore the error in the cumulant Cm grows 



with (n) 



(m+l)/2 



Overall, a reasonable agreement be- 



tween theory and experiment is obtained. Universal os- 
cillations of the cumulants^^ highlight the difference be- 
tween the distribution and a Gaussian distribution for 
which Cn = for n > 2 and provide additional informa- 
tion about the probability distribution. An interesting 
piece of information is the position of the zeros of the 
generating function (ZGF) in the complex plane which, 
according to Abanov et al.^, is expected to be on the 
negative real axis for non-interacting systems. This is 
interesting as the strong interactions in p-type QDs may 
cause deviations from single-particle physics. However, 
it is difficult to extract any useful information directly 
from normal cumulants as the poles of the cumulant gen- 
erating function are displaced from the real axis by con- 



We have demonstrated time-averaged as well as time- 
resolved charge detection in a p-type GaAs QD. The ex- 
tracted tunnelling rates suggest the presence of a dense 
spectrum of excited states in the dot contributing to 
transport. The full counting statistics of the QD is stud- 
ied and shown to follow the two-state Markovian BN 
model in- spite of multi- level transport. This result and 
the absence of QPC back-action are interpreted in terms 
of a strong energy relaxation of holes in the QD which 
also ensures uni-directional transport within the band- 
width of our measurement. 



Appendix 

Hidden states are a set of internal states of a system 
that are indistinguishable from the perspective of the 
charge detector. These could be some two-level systems 
in the barriers which affect the rates, some excited energy 
states within the dot or even some internal states of the 
charge detector. In this Appendix we consider a sim- 
ple model of the presence of an excited state in addition 
to the ground state of the dot, showing that when the 
hidden state is traced out, the dot appears to be obey- 
ing non-Markovian statistics. However, in the presence 
of strong relaxation, Markovian statistics is recovered. 
Similar problems have been considered by Belzi^^ and 
Flindt et al.^^ The starting point is the master equation 
for the three state system 
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where F-^ = Ff^ ± Ff^ and a similar expression for P^^^ . 
R^^ is the relaxation rate, E^^ is the rate of photon- 
assisted excitations and z is the complex-value counting 



7 



THHII (b) 150 




.... .... .^""iS^^J 


'A — ^ 


-e-C8/ci 
-e-C9/ci 
-e-cio/ci 

C11/C1 

-e-ci2/ci 


:„ llilllll 



10 10 10 

Mean Event Number 




10 10 
Mean Event Number 



-e- positive 


F12 


-e- negative 


^(^^ F11 




^^^^^*^* F10 




F9 















10 10 
Mean Event Number 



10 10 
Mean Event Number 



FIG. 5: (a)-(b) First twelve normalized cumulants of charge 
transfer as a function of time (mean number of events in time 
traces with a given length cut from a very long time trace) 
calculated from the data collected at point A in Fig.|2ja) along 
with the predictions of BN model which fits remarkably well 
to the data. The amount of statistics decreases for longer 
time traces causing the size of error bar to increase, (c)- 
(d) First twelve factorial cumulants (FCs)^^ calculated from 
the normal cumulants. FCs of different order alternate sign 
and grow exponentially with time and therefore it is more 
convenient to follow their trend in logarithmic scale as shown 
in (e)-(f). No zero-crossing oscillations in FCs are observed 
consistent with a two-level Markovian model. The white 
arrows point to the development of faint non-zero-crossing 
oscillations appearing in higher order FCs probably due to 
finite statistics. 



fielcP^ = Pi{z^t) for i = O^g^e and standard results 
are recovered for 2; = 1). The charge detector is insensi- 
tive to the particular state of the dot, therefore writing 
Pi = Pg Pe^ we have 



Po ^ 






out 


^} 1 
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^ out 
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. (3) 



Concentrating on the visible subspace by defining 



w : 



Mv 



M 



-r+ 

m 

r.+„ 



^-'^out 
^ out 



(4) 



where v = (po Pi)" , it can be seen that deviations from 
Markovian statistics (w = 0) are caused by the popula- 
tion of the excited state 



W : 



Pe- 



(5) 



This deviation is also proportional to F^^^ and it vanishes 
for Fq^^ = Fq^^ as the two states become statistically 
indistinguishable. In the limit of i?^^ > rf^,£;^^, this 
population vanishes {pe 0) in the steady-state, and 
the Markovian solution is recovered. This can be seen by 
using Eq.[2]to eliminate Pe from the previous equation 



w + (F^,,+i?^^)w = F- 



^ m 

-F? 



(6) 



In the limit i?^^ > Tf^.E'^^, the right hand side can be 
neglected and w = solves the resulting equation. 
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